Universal elastic and plastic effects in the particle caging of polymers and glassforming liquids 
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The peculiar intermittency of the single-particle motion close to the glass transition follows from the contin- 
uous (elastic) trapping into the neighbors cage and the subsequent (plastic) escape. The universal correlation 
of the residual elasticity and the incipient plasticity on the picosecond time scale with the cage restructuring is 
revealed by extensive simulations and validated by comparison with experimental data spanning about eighteen 
decades in relaxation time and a wide range of fragilities. Universality is lost on sub-picosecond time scales 
where the elastic energy due to non-affine microscopic displacements is still stored. 

PACS numbers: 64.70.P-, 62.20.de, 66.20.-d 
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If liquids avoid crystallization during cooling or compres- 
sion, their viscosity grows by orders of magnitude and they fi- 
nally freeze into a microscopically disordered solid-like state, 
a glass yj]. On approaching the glass transition (GT), par- 
ticles tend to be trapped in transient cages formed by their 
nearest neighbours. The particles rattle about in the cage on 
picosecond time scales and are later released with average es- 
cape time, or structural relaxation time r Q , increasing from a 
few picoseconds to thousands of seconds through the transi- 
tion JIII2]. After the release, particles undergo short periods of 
larger displacement terminated by a new trapping [3|,|4J]. Ow- 
ing to this stop-and-go process, the single-particle dynamics 
becomes strongly intermittent |5J]. 

The fast rattling motion of particles during the trapping 
periods has strong analogies with the oscillatory elastic be- 
haviour of particles in crystalline and amorphous solids. Sev- 
eral "elastic" models of the structural relaxation and the vis- 
cous flow have been developed since long time ago [6|] and 
over the years (for excellent reviews see 1 71 l al )- See also re- 
cent experimental SOU], theoretical iflUfL?!! and computa- 
tional | T(], IH [l4ll work in a wide class of materials includ- 
ing colloids, liquids, polymers and metallic alloys. The key 
prediction of the elastic models is that the viscosity rj and 
the structural relaxation are activated with energy cost, due 
to the cage breaking, being proportional to the instantaneous, 
or infinite-frequency, shear modulus Goo QUI]: 

A K 

log 77/770, logT a /r = — ; AE = G X V* (1) 
k B T 

where V* is an activation volume with negligible temperature 
dependence. EqQ] relies on two underlying assumptions: i) 
only the trapping period, and not the subsequent escape pro- 
cess, is pertinent to relaxation, ii) continuum elasticity theory 
applies to microscopic length scales as well. Both assump- 
tions are debatable in that: i) departing from GT, trapping 
is less effective and the role of the escape process increases, 
ii) after a homogeneous large-scale deformation, the discrete- 
ness of the matter leads to inhomogeneous, non-affine fast dis- 
placements of the particles within the cage 1 15 , JiJ with the 



subsequent drop of the modulus from Goo to the finite fre- 
quency plateau (or "relaxed") modulus G r II 1 711 . Virtually, all 



the experiments validating Eq[T]actually measure G r since the 
elastic response is probed at frequencies much lower than the 
THz range where the crossover from G r to Goo is located 
[ lT|. Then, the approximation Goo — G r is usually assumed 
17-flol 121. Surprisingly, in spite of the delicate conjectures 
underlying the popular elastic models a thorough test of Eq[T] 
is missing. 

The present Letter reports striking evidence by extensive 
molecular-dynamics (MD) simulations that G r , but not Goo, 
exhibits strong universal correlations with the short-time (pi- 
cosecond) mean square displacement (ST-MSD) of the kinetic 
unit (u 2 ). The latter quantity, which will be related to specific 
details of the potential-energy landscape (PEL) 1181 . reveals a 
simple picture of the intermittent dynamics expressed by the 
quasi-elastic relation (in reduced MD units): 
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It is shown that the first term on the rhs is the inverse ST- 
MSD during the trapping period, whereas (u^) is associated 
to the incipient escape from the cage. Eq|2]improves the usual 
elastic approximation with (m^,)^ 1 = and G r = Goc [7 



11 



1911 . The combination of Eq|2]with previous results [1^ 
predicts the scaling of r a , or 77, and G r /T. It is shown that 
the resulting master curve fits the experimental results over 
about eighteen decades in relaxation time (—12 < log r a < 6) 
and a wide range of fragilities (20 < m < 115) with one 
adjustable parameter, to be compared with the reported range 
of validity of Eq. [U(— 6 < log r Q < 2 (|7|]) with two adjustable 
parameters (V* and To or 770) and the identification Goc = G r . 

We performed extensive molecular-dynamics (MD) simula- 
tions of a melt of fully-flexible linear chains of M soft spheres 
(monomers, N ~ 2000 in total) by changing the tempera- 
ture T, the number density p, the chain length and the in- 
teracting potential between non-bonded monomers. Details 
about the model are given elsewhere II1411 and summarized in 
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Supplemental Material (SM) |20ll where all the investigated 
states (~ 100) are also listed. MD results are in reduced units 
(Boltzmann constant kg = 1). 

We start by comparing both the single-particle and the col- 
lective dynamics. The former is characterized via the in- 
coherent intermediate scattering function (ISF) F s (q,t) — 
N ~ 1 (I2f=i exp{-iq ■ [rj(t) - r,(0)]}) at q rnax , the q-vector 
of the maximum of the static structure factor ||2|] . The relation 
Fs (qmax ,T a ) = 1/e defines the structural relaxation time t q . 
The collective dynamics is described by the transient elastic 
modulus of a volume V, G(t), being expressed by the corre- 
lation function: 



G(t) 



V 



{Vxy(t )a xy (t + t)). 



a xy is the off-diagonal component of the stress tensor : 



iff 
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(3) 



(4) 



where v a k, F a ki, r a ki are the a components of the velocity of 
the fcth monomer with mass m, the force between the fcth and 
the Ith monomer and their separation, respectively. Note that 
G(0) = Goo & 

FigQ]shows the characteristic two-step decay of both G{t) 
(top) and ISF (bottom). The drop of both G(t) and ISF over 
times t ~ 1, corresponding to few picoseconds in real time 
JI4I1 . follows from the exploration of the cage by the trapped 
monomer. This results in the fast, large drop of G(t) from Goo 
to G r , i.e. the plateau height measured by most viscoelastic 
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experiments 111 111 . It seems proper to deepen the discussion 
about this point. In crystals with one particle per unit cell un- 
der athermal conditions the equality G r = Goo holds IU5I1 . In 
such a case, the total force acting on each particle, F tot van- 
ishes during the shear strain. Differently, even under ather- 
mal conditions, in crystals with multi-atom unit cell fl 1 5fl and 
amorphous systems [16] the total force on each particle un- 
balances (Ftot 7^ 0) after the fast, homogeneous, affine shear 
strain. At finite temperature the mechanical equilibrium is re- 
stored rapidly by the non-affine particle displacements within 
their cages, leading to considerable loss of elastic energy 11611 - 
Note that the differences between G r and Goo are reduced but 
not suppressed in the glass for T < T g since local cage rat- 
tling is still present 12211 . 

The fact that G(t) vanishes at t ~ r Q (FigQ]top) shows that 
the plateau G r is controlled by the local structure lfl7ll . To 
provide more insight, the configurations explored by the sys- 
tem via the true dynamics are sampled in time and "quenched" 
to the so-called inherent structures (IS) by locally minimizing 
the potential energy in configurational space, i.e. enforcing 
mechanical equilibrium (F tot — 0) [18]. The IS time series 
defines the "inherent" dynamics of the system. FigQ] shows 
that the true and the inherent decay of G(t) differ at very 
short times but coincide for times longer than t ~ 1. This 
confirms that the initial drop of G(t) follows from mechani- 
cal equilibration. The insets of FigfUtop) provide information 




1 

0,8 



a 0,6 



3 



0,4 
0,2 







-*>J 1 1 1 1 1 1 1 1 


1 1 1 1 1 


T \ \ 




— 0.3 \ \ 




— 0.34 \ 




— 0.4 \ 

V- 1 





-2 -1 



1 2 3 

log(t) 



FIG. 1: (Color online) MD results of the relaxation properties of 
linear chains with M = 3, p = 0.984 at the indicated temperatures. 
Non-bonded monomers interact with Lennard-Jones potential 12011 . 
Top Panel: transient elastic modulus G(t). Note that G(0) = Goo- 
The relaxed shear modulus is defined by G r = G(t*), t* is defined 
in the text. Insets show the true and the quenched (inherent) radial 
pair distribution functions in the regions of the first and the second 
neighbor shells at T = 0.3. The peak at r ~ 0.97 corresponds to 
the adjacent bonded monomers. The fine structure of the peak of 
the second neighbor shell (1.5 < r < 2.5) is discussed in the text. 
Bottom Panel: intermediate scattering function F s (q max , t). Dots 
mark the structural relaxation time r Q . 



on the IS structures by comparing the actual and the inherent 
radial distribution functions. The latter, having removed the 
vibrational cage rattling, is much more detailed. In particular 
the peak corresponding to the second neighbor shell is split in 
three sub-peaks at r ~ 1.89,2.05,2.19. The side sub-peaks 
are well-known signatures of local icosahedral order in atomic 
Lennard-Jones liquids at r ~ 1.89, 2.15 12311 . The central sub- 
peak is characteristic of the present molecular system and is 
assigned to linear trimers of two bonded and one non-bonded 
monomers (not shown). 

The first escapes of the monomer from the cage occur at 
very short times. By observing the monomer MSD (r 2 (t)), 
it was noticed that the quantity 9 log (r 2 (t)) / 'dlogt shows a 
well-defined minimum at t* = 1.023 signaling early detrap- 
ping (i* is independent of the physical state in the present 
model) J3. We define (u 2 ) = (r 2 (t*)) and G r = G(t*) 
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high correlations of G r with both (u 2 ) and r Q are observed 
(FigfJJ?). In particular, the inset shows the remarkable linear 
scaling law, Eqf2] |25]]. The constancy of the parameters of 
Eq. |2] V* and (u 2 ^), is ascribed to the limited changes of 
the local structure due to the high packing of the investigated 



states 124 12611 . The combination of Eqf2]and Eq [5] yields a 
quadratic master curve between logr a and G r /T which well 
agrees with the MD results (Fig[2t>) and is not recovered by re- 
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140 160 placing Goo with G r in Eqfl] FigFJ]} also shows that (u 2 ^)' 



cannot be neglected in Eq|2] i.e. far from GT the elastic re- 
sponse competes with, i.e. is often interrupted by, the escape 
process. However, the elastic term dominates close to GT, 
thus recovering the usual elastic approximation J3, EI , 19 ] . 

To provide insight into Eq. [2] we resort to the PEL view- 
point, where IS are grouped to form distinct metabasins (MB) 
separated by high energy barriers L 1 Si] - In this context one 
sees that the elasticity due to trapping is set by intra-MB dy- 
namics, whereas the trapping release is accounted for by the 
inter-MB dynamics. More quantitatively, we consider the re- 
turn distance (squared) (Sr 2 ) = N" 1 X)ili(l r «(^) — r t\ 2 )^ 
involving the separation between the instantaneous and the 
quenched IS positions, thus providing an estimate of the aper- 
ture of that portion of the basin encountered by the trajec- 
tories [18]. Remarkably, Sr 2 , differently from the instanta- 
neous displacement involved in the definition of (u 2 ), is a 
stationary random variable and then (Sr 2 ) is a characteris- 
tic constant quantity of each physical state. It is found that 
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FIG. 2: (Color online) Panel (a): The structural relaxation time r Q 
and the ST-MSD (u 2 ) (inset) versus the ratio G^jT from MD sim- 
ulations. Panel (b): Same as in panel (a) with the reduced elas- 
tic modulus G r - The open triangles and the dashed heavy line in 
the inset are the inverse return distance (Sr 2 )" 1 of the states in 
Fig[]](plus T = 0.6,1) QJ] and Eq|2]with V* = 0.77(2) and 
(u 2 )r X1 = 0.192(4), respectively. The dashed heavy and thin lines 
in the main panel are the combination of Eq. [5] with Eqf2] by ei- 
ther including or not the term (u 2 )^ of Eqfl] respectively. Note 
that replacing Goo by G r in Eq[T]leads to a linear master curve dif- 
fering from the MD results. The investigated states, obtained by 
changing T, p, M and the interacting potential between non-bonded 
monomers, are listed in SM [ 20] . 



as the ST-MSD [14] and the residual elastic modulus (see 
FigfUtop)), respectively. It is known that strong correlation 
exists between the structural relaxation and (u 2 ) [24], result- 



(5) 



ing in the universal master curve [ 1 

logr Q = a + fiiu 2 )- 1 +~i{u 2 Y 



with a = -0.424(1), /3 = 2.7(1) • 1(T 2 , 7 = 3.41(3) • 1(T 3 
in MD units. 

We are now in a position to investigate the correlation of 
the elastic properties with (u 2 ) and then t q . Fig f2] compares 
the latter quantities with both Goo (top) and G r (bottom). It 
is quite apparent that Goo, the central quantity of the standard 
elastic models ||7t|8|], poorly correlates with both the structural 
relaxation and even (u 2 ) despite the latter is evaluated at very 
short times, few picoseconds in real time. Differently, very 



{Sr 2 ) = 0.625(u 2 ) (FigfJl', inset), i.e. (u 2 ) is an excellent 
surrogate of the return distance (the proportionality constant 
is close to 0.5, expected for uncorrelated deviations of the in- 
stantaneous positions from the IS). This identification allows 
us to interpret the two terms of Eq. (u 2 ^) and the "elastic" 
ST-MSD k B T/G r V* dominating close to GT, as the average 
MB full width (or approximately the MSD to reach a different 
MB), and the typical MSD within a MB, respectively. As an 
example, for a binary mixture the typical MSD within a MB 
(~ 6 ■ 10 -2 ) and to reach a different MB (~ 0.2) are rather 
comparable, at the same Tq,(~ 4 • 10 3 ), to the "elastic" ST- 
MSD (~ 4.5 • 10~ 2 ), and (u 2 )oc ^ 0.192 of the monomers, 
respectively. 

Are the MD results robust enough to be compared with the 
experiments ? The fact that Eqf5] assessed in the weakly su- 
percooled regime attainable by MD simulations, holds true 
down to GT (logT QS = 13.5 ± 0.5, (u 2 g ) = 0.0166) evi- 
dences that the relevant information on the coupling between 
the intermittent motion and the long-time relaxation is present 
at high temperatures [14]. This motivates us to extend the 
validity of Eqf2]down to GT and get the ratio G rg /T g since 
(u 2 g ) = 0.0166 [14]. By plugging Eq. fjinto Eq. Hand defin- 
ing the reduced variable X = G r T g /(G rg T), one obtains a 
quadratic form as candidate for the universal master curve: 



logr Q = 6 + eX + (j)X 2 



(6) 



with S = -11.70(1), e = 3.4(2), <j> = 
to get the familiar log r a = 2 at GT). 



10.3(8) (6 is set so as 
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FIG. 3: (Color online) Scaling of the structural relaxation time and 
viscosity in terms of the reduced variable X = G r T g / (G rg T). 
The continuous line is the quadratic form Eq(6] The dashed 
straight line is the linear form EqQ] having replaced Goo with 
G r - The numbers in parentheses denote the fragility of the sys- 
tems (data sources in SM [20]). Data are considered only if the 
corresponding elastic moduli are measured at (angular) frequency 
u > max {0.1/T a , 2ty ■ 1 MHz }. Inset: \ 2 of the best-fit of the 
experimental data with both the quadratic form Eq[6] (allowing only 
the vertical shift of the experimental data [20]) and the modified EqQ] 
(two adjustable parameters: the slope and the vertical shift). The fit 
spans the variable range from log t^ /to up log r« /to = 6. 



Fig|3] shows the comparison of Eq|6]with several systems 
spanning a wide range of fragilities and elastic moduli mea- 
sured at sufficiently high frequencies to allow comparison 
with t q or -q over wider ranges (elastic response is observed 
if u)T a ^> 1) . The experimental data are vertically shifted 
to adjust the conversion factor between MD and real units 
and ensure best-fit with Eq|6] No other adjustment is done. 
Data sources and related shifts are listed in SM [20]. Figj3] 
shows that the modified EqQ] with Goo being replaced by 
G r results in a poorer fit in spite of the additional adjustable 
parameter (the slope). To get further insight, we calculated 
the chi-square of the two fits over the variable range from 
logriT /tq up log Ta™ ^ /t = 6 (Fig |3] inset). It is seen 
that the two fits are equivalent in the range —4 < log r Q < 6 
but Eq|6]performs much better in the range — 12 < log r a < 
—4. This is expected in that our analysis accounts for the 
competing effects of elastic trapping and anelastic escape, a 
feature which is missing in standard elastic models J3,lltl- 

In conclusion, we showed that plasticity cannot be ne- 
glected even on picosecond time scales where elastic effects 
seem to dominate the intermittent dynamics. The emerg- 
ing quasi-elastic picture reveals the universal correlations be- 
tween the residual elasticity left by intra-cage relaxation and 
the long-time relaxation due to cage breaking in a range of 
about eighteen decades in relaxation times, as well as clear 
connection to the underlying PEL. 
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